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Abstract 

We discuss the physics interest and the experimental feasibiUty for detecting 
high-mass dimuon pairs using the planned 50-GeV Proton Synchrotron (PS) at the 
KEK/JHF and JAERI/NSP joint accelerator project. The Drell-Yan measurement 
oi p + d versus p + p aX 50 GeV will provide unique information on the flavor asym- 
metry of proton's up and down sea-quark distributions in the large-x region. A 
study of the nuclear dependences of Drell-Yan cross sections can reveal the modi- 
fication of antiquark distributions in nuclei. Furthermore, the effect of energy loss 
for fast partons traversing nuclear medium could also be sensitively measured. If 
polarized proton beam becomes available at the 50-GeV PS, unique information on 
the sea-quark polarization could be obtained. Study of heavy quarkonium produc- 
tion at the 50-GeV PS can set important constraints on the mechanism of vector 
meson productions. Using a prototype dimuon spectrometer, we have simulated the 
sensitivities for a variety of measurements. 
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1 Introduction 



One of the most active areas of research in nuclear and particle physics during the last 
several decades is the study of quark and gluon distributions in the nucleons and nuclei. 
Several major surprises were discovered in Deep-Inelastic Scattering (DIS) experiments 
which profoundly changed our views of the partonic substructure of hadrons. In the early 
1980's, the famous 'EMC effect found in muon DIS provided the first unambiguous evi- 
dence that the quark distributions in nuclei are significantly different from those in free 
nucleons 0. More recently, surprising results on the spin and flavor structures of the 
nucleons were discovered in DIS experiments. The so-called 'spin crisis', revealed by the 
disagreement between the prediction of the Ellis- Jaffe sum rule and the polarized DIS 
experiments, has led to extensive theoretical and experimental efforts to understand the 
partonic content of proton's spin . Subsequently, the observation |Q of the violation of 
the Gottfried sum rule in DIS revealed a surprisingly large asymmetry between the up 
and down antiquark distributions in the nucleon, shedding new light on the origins of the 
nucleon sea. 



The partonic structure of nucleons and nuclei can also be measured with hadronic 
probes. A powerful tool for such studies is the Drell-Yan process 0, in which a quark 
annihilates with an antiquark forming a virtual photon which subsequently decays into a 
lepton pair. The proton-induced Drell-Yan process is of particular interest, since it can 
be used to extract antiquark distributions of the target nucleon and nuclei. This provides 
information complementary to what can be obtained in DIS, which is sensitive to the sum 
of the quark and antiquark distributions. 



The usefulness of the Drell-Yan process as a tool for probing antiquark distributions 
has been well demonstrated by a series of Fermilab dimuon production experiments . In 
particular, the Drell-Yan cross section ratios for p + d versus p + p led to a direct measure- 
ment of the d/u asymmetry as a function of Bjorken-x [^, |^. Furthermore, the nuclear 
dependence of the Drell-Yan cross sections showed no evidence for antiquark enhancement 
in heavy nuclei |rD|, in striking disagreements with predictions of some theoretical models 
which were capable of explaining the EMC effect. 



The 50-GeV Proton Synchrotron (PS) offers a unique opportunity to extend existing 
measurements of antiquark distributions to much larger values of Bjorken-x. Such infor- 
mation is crucial for understanding the origins of flavor asymmetry in the nucleon sea, and 
for illuminating the nuclear environment effects on parton distributions. Moreover, Drell- 
Yan measurements using polarized proton beam on polarized target at the 50-GeV PS will 
provide a first determination of the spin-dependent antiquark distribution at an x region 
not accessible in the RHIC-spin program. Indeed, the flavor asymmetry of polarized sea 
quark distributions, predicted to be very large by certain theoretical models |n], |12|, can 
be directly measured. 
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A detailed study of the nuclear dependence of Drell-Yan cross sections at 50 GeV 
could also lead to a first observation of the coherent partonic energy-loss effects predicted 
recently by Baier, Dokshitser, Mueller, Peigne, Schiff (BDMPS) [1^ and by Zakharov ||T^ . 
These authors studied the radiative energy loss (through gluon emission) of high energy 
partons passing through hot and cold hadronic matter. The partonic energy-loss effect is 
the QCD analog of the Landau-Pomeranchuk-Migdal (LPM) QED effect |jl5|, |16[ predicted 
over 40 years ago and confirmed only recently at SLAG fl^. A number of surprising ef- 
fects were obtained by BDMPS and Zakharov. First, the partonic energy loss in a hot 
QCD plasma is predicted to be much larger than in cold matter. This suggests that an 
anomalously large energy loss of jets produced in relativistic heavy-ion collisions could be 
a signature for Quark-Gluon Plasma formation. Second, the radiative energy loss is pre- 
dicted to be proportional to L^, where L is the path length of hot or cold nuclear matter 
traversed by the partons. This curious result is contrary to the conventional wisdom that 
energy loss depends linearly on L, and it reflects the quantum-mechanical interference 
effect from several contributing diagrams. 



An attempt to search for partonic energy-loss effects in cold matter was made recently 



via the study of nuclear dependence of Drell-Yan cross sections at 800 GeV |I8[ . Only an 
upper limit for partonic energy loss was determined. A much more sensitive study can be 
made at lower beam energies, where the fractional energy loss /S.E / E will be larger. At 
the 50-GeV PS, the effect is expected to be much enhanced and indeed one could even 
examine whether the nuclear effect follows the L or dependence. 



While logarithmic scaling violation is well established in DIS experiments, no clear 
evidence for scaling violation has been seen in Drell-Yan process. The 50-GeV PS pro- 
vides an interesting opportunity for unambiguously establishing scaling violation in the 
Drell-Yan process |]l9l. For given values of xi and X2 (Bjorken-x for the projectile and 
target partons, respectively), scaling- violation is expected to cause roughly a factor of two 
increase in the Drell-Yan cross sections when proton beam energy is decreased from 800 
GeV to 50 GeV. It appears quite feasible to establish scaling violation in the Drell-Yan 
process with future data from the 50-GeV PS. 



Detection of high- mass dileptons at the 50-GeV PS will also allow a study of J/\E' and 
^' production. Existing data on proton-induced charmonium production are mostly lim- 
ited to the energy range 150 GeV < Ep < 800 GeV. A comparison of 50 GeV charmonium 
production data with existing data will further improve our knowledge on the produc- 
tion and propagation of charmonium in the nuclear medium. Many different effects which 
could affect the production of charmonium in nuclear medium, such as nuclear shadowing, 
partonic energy loss, final-state interaction with comoving gluons or hadrons, will have 
different beam energy dependences 21|. Therefore, a systematic study of charmonium 



production in p — p,p — A and A — A collisions at 50 GeV would be extremely valuable 
for disentangling various effects. Only after the mechanisms for charmonium production 
are well understood could J/\Ef-suppression be used confidently as a signature for Quark- 
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Gluon Plasma formation in relativistic heavy-ion collisions 



Many of the proposed studies for Drell-Yan process at the 50-GeV PS could also ben- 
efit from the J/\E' and \E'' data. Unlike the Drell-Yan which is an electromagnetic process, 
quarkonium production is a strong-interaction process involving gluon-gluon fusion and 
quark-antiqaurk fusion. Comparison between the Drell-Yan and quarkonium production 
data will further elucidate various aspects of parton distributions in nucleons and nuclei, 
and of the propagation of partons in nuclei. As an example, charmonium production with 
polarized proton beam at the 50-GeV PS might provide interesting information on the 
gluon distributions at large Bjorken-x, which is essentially unknown. 

In this paper, we discuss the physics interest and the feasibility for making precise 
measurements of high-mass dimuons at the 50-GeV PS. In Section 2 the physics motivation 
for various measurements will be discussed in some details. A preliminary design study 
of a dimuon spectrometer for the 50-GeV PS will be presented in Section 3. A summary 
is given in Section 4. 



2 Physics Issues 

2.1 Overview of high-mass dilepton production 

Detection of high-mass dileptons produced in high-energy hadronic interactions has a 
long and glorious history. The charm and beauty quarks were discovered in the 1970's 
via the dilepton decay modes of J/^ and T resonances. These quarkonium states are 
superimposed on a dilepton continuum known as the Drell-Yan process M. The Drell- 



Yan data has been a source of information for the antiquark structure of the nucleon ||24 . 
Furthermore, Drell-Yan production with pion and kaon beams has yielded the parton dis- 
tributions of these unstable particles for the first time. A generalized Drell-Yan process 
was also responsible for the discovery of the W and Z gauge bosons in the 1980 's. 

To lowest order, the Drell-Yan process depends on the product of quark and antiquark 
distributions in the beam and target as 

, , = 7{ $]e^[g„(xi)ga(x2) + ga(a;i)ga(a;2)]- (1) 

dXidX2 ys X1X2 a 

Here qa{x) are the quark or antiquark structure functions of the two colliding hadrons 
evaluated at momentum fractions Xi and X2- The sum is over quark flavors, and s is the 
center-of-mass energy squared. 

The kinematics of the virtual photon — longitudinal center-of-mass momentum 
transverse momentum pj- and mass — are determined by measuring the two-muon 
decay of the virtual photon. These quantities determine the momentum fractions of the 
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two quarks: 

Xp = vl I vl^"^""^ = Xi — X2 (2) 

M2 = X1X2S (3) 

where pj is the virtual photon center-of-mass longitudinal momentum and pj'"^""^ is the 
maximum value it can have. 



Although the simple parton model enjoyed considerable success in explaining many 
features of the early data, it was soon realized that QCD corrections to the parton model 
were required. The inclusion of the NLO diagrams for the Drell-Yan process brings excel- 
lent agreement between the calculations and the data. As an example, Figure |I| shows the 
NA3 data [||] at 400 GeV, together with the E605 and E772 |3 data at 800 GeV. 
The solid curves in Figure |1| correspond to NLO calculation for 800 GeV p + d {^/s = 38.9 
GeV) and they describe the NA3/E605/E772 data well. This shows that the mechanism 
for Drell-Yan process is well understood theoretically, and quantitative information on 
the parton distributions can be reliably extracted via this process. 



To gain sensitivity to the antiquark distribution of the target, one chooses a proton 
beam and selects the kinematic region of positive xp and large xi. In this limit the 
contribution from the second term in Eq. |^ is small and the first term is dominated by the 
u{xi) distribution of the proton. Under these circumstances, the ratio of the cross sections 
for two different targets, X and Y, which have Ax and Ay nucleons is approximately the 
ratio of the u{x2) distributions: 

1^ {dxidx2) ^ U^{X2) 
Ay \dx\dx2 J 

In this relation the cross sections are defined per nucleus but the parton distributions are 
conventionally defined per nucleon. 



(4) 



Eq. ^ demonstrates the power of Drell-Yan experiments in determining relative an- 
tiquark distributions. This feature was explored by recent Fermilab experiments using 
800 GeV proton beams 0. The 50-GeV PS provides a unique opportunity for extending 
the Fermilab measurements to larger X2 {x2 > 0.25). For a given value of xi and X2, the 
Drell-Yan cross section is proportional to 1/s (see Eq. |l]). Hence the cross section at 50 
GeV is roughly 16 times greater than that at 800 GeV (The price one pays at lower beam 
energies is that one has limited reach for small X2, which could best be studied at higher 
energies). Furthermore, to the extent that the radiation dose scales as beam power, one 
can take ~ 16 times higher beam flux at 50 GeV relative to 800 GeV. The combination 
of these two effects could lead to two orders of magnitude improvement in the statistics 
at high X2 over previous Fermilab experiments. 
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2.2 Scaling violation in Drell-Yan process 



The right-hand side of Eq. |I| is only a function of xi,X2 and is independent of the beam 
energy. This scahng property no longer holds when QCD corrections to the Drell-Yan 
process are taken into account. While scaling violation is well established in DIS exper- 
iments, it is not confirmed in Drell-Yan experiments at all. No convincing evidence for 



scaling violation is seen As discussed in a recent review [0], there are mainly two 
reasons for this. First, unlike the DIS, the Drell-Yan cross section is a convolution of 
two structure functions. For proton-induced Drell-Yan, one often involves a beam quark 
with xi > 0.1 and a target antiquark with X2 < 0.1. Scaling violation implies that the 
structure functions rise for x < 0.1 and drop for x > 0.1 as increases. Hence the 
effects of scaling violation are partially cancelled. Second, unlike the DIS, the Drell-Yan 
experiment can only probe relatively large Q^, namely, > 16 GeV^ for a mass cut of 4 
GeV. This makes it more difficult to observe the logarithmic variation of the structure 
functions in Drell-Yan experiments. 

The 50-GeV PS provides an interesting opportunity for unambiguously establishing 
scaling violation in the Drell-Yan process. Figure ^ shows the predictions for p + d at 
50 GeV. Scaling violation causes a factor of two increase in the Drell-Yan cross sections 
when the beam energy is decreased from 800 GeV to 50 GeV. It appears quite feasible 
to establish scaling violation in Drell-Yan with future dilepton production experiments at 
the 50-GeV PS. 



2.3 d/u asymmetry of the proton 



From neutrino-induced DIS experiments, it is known that the strange quark sea in the 
nucleon is roughly a factor of two less than the up or down quark sea The lack 

of SU(3) flavor symmetry in the nucleon sea is attributed to the much heavier mass of 
the strange quark. Until recently, it had been assumed that the distributions of u and d 
quarks were identical. Although the equality of u and d in the proton is not required by 
any known symmetry, this is a plausible assumption for sea quarks generated by gluon 
splitting. Because the masses of the up and down quarks are small compared to the con- 
finement scale, nearly equal numbers of up and down sea quarks should result. 



The assumption of u{x) = d{x) can be tested by measurements of the Gottfried 
integral 0, defined as 



ix,Q' 



^2"(x,Q2)1 /:j;dx = \ + l 

J 6 6 JO 



dp{x) 



dx. 



(5) 



where and are the proton and neutron structure functions measured in DIS ex- 
periments. Under the assumption of a symmetric sea, u = d, the Gottfried Sum Rule 
(GSR) 1^, Ig = 1/3, is obtained. The most accurate test of the GSR was reported in 
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1991 by the New Muon Collaboration (NMC) 0], which measured and F2 over the 
region 0.004 < x < 0.8. They determined the Gottfried integral to be 0.235 ± 0.026, sig- 
nificantly below 1/3. This surprising result has generated much interest, and it strongly 
suggested that the assumption u = d should be abandoned. Specifically, the NMC result 
implies 

/ rfx = 0.148 ±0.039. (6) 

JO L J 

Eq. ^ shows that only the integral oi d — u was deduced from the DIS measurements. The 
X dependence oi d — u remained unspecified. 

The proton-induced Drell-Yan process provides an independent means to probe the fla- 
vor asymmetry of the nucleon sea ||32| . An important advantage of the Drell-Yan process 



is that the x dependence of d/u can be determined. The Fermilab E772 collaboration 
compared the Drell-Yan yields from isoscalar targets with that from a neutron-rich (tung- 
sten) target, and constraints on the nonequality of u and d in the range 0.04 < x < 0.27 



were set. More recently, the CERN experiment NA51 |3J| carried out a comparison of the 
Drell-Yan muon pair yield from hydrogen and deuterium using a 450 GeV/c proton beam. 
They found that u/d = 0.51 ± 0.04 ± 0.05 at (x) = 0.18, a surprisingly large difference 
between the u and d. 

A Drell-Yan experiment (E866), aiming at higher statistical accuracy and wider kine- 
matic coverage than NA51, was recently completed p, ^ at Fermilab. This experiment 
also measured the Drell-Yan muon pairs from 800-GeV/c protons interacting with liquid 
deuterium and hydrogen targets. Eq. |^ shows that the Drell-Yan cross section ratio at 
large xp is approximately given as 

apYip + d) _ 1 A , d{x2) \ , . 

2aDY{p + p) ^ 2 1^ u{x2)J' ^ ' 

Values for d/u were extracted by the E866 collaboration at = 54 GeV^/c^ over the 
region 0.02 < x < 0.345. These are shown in Figure ^ along with the NA51 measure- 
ment. For X < 0.15, d/u increases linearly with x and is in good agreement with the 
CTEQ4M ||28[ and MRS(R2) |3^ parameterizations. However, a distinct feature of the 



data, not seen in either parameterization, is the rapid decrease toward unity oi d/u be- 
yond X = 0.2. 

The d/u ratio, along with the CTEQ4M values for d + u, was used to obtain d — u 
(Figure ||). Being a flavor nonsinglet quantity, d{x) — u{x) is decoupled from gluon dis- 
tribution. Since perturbative processes have negligible contribution to d/u asymmetry, 
d{x) — u{x) essentially isolates the contribution from non-perturbative effects. From the 
results shown in Figure one can obtain an independent determination of the integral 
of Eq. § E866 finds 0.100 ± 0.007 ± 0.017, consistent with, but roughly 2/3 of the value 
deduced by NMC. 
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As early as 1983, Thomas pointed out that the virtual pions that dress the pro- 
ton will lead to an enhancement of d relative to u via the (nonperturbative) "Sullivan 
process." Sullivan previously showed that in DIS virtual mesons scale in the Bjorken 
limit and contribute to the nucleon structure function. Following the publication of the 
NMC result, many papers treated virtual mesons as the origin of the d/u asymmetry 



for a recent review). Here the TT'^{du) cloud, dominant in the process p —>■ 



see 



leads to an excess of d sea. 



n, 



A different approach for including the effects of virtual mesons has been presented by 
Eichten et al. and further investigated by other authors |^ , In chiral perturba- 
tion theory, the relevant degrees of freedom are constituent quarks, gluons, and Goldstone 
bosons. In this model, a portion of the sea comes from the couplings of Goldstone bosons 
to the constituent quarks, such as m ^ dn'^ and d uti~ . The excess of d over u is then 
simply due to the additional valence u quark in the proton. 

The X dependence oi d — u and d/u obtained by E866 provides important constraints 
for theoretical models. Figure ^ compares d{x) — u{x) from E866 with a virtual-pion 



model calculation, following the procedure detailed by Kumano . A dipole form, with 
A = 1.0 GeV for the ttNN form factor and A = 0.8 GeV for the vrA^A form factor, was 
used. A is the cutoff parameter for the pion form factor. Figure § also shows the predicted 
d — u from the chiral model The chiral model places more strength at low x than 
does the virtual-pion model. This difference reflects the fact that the pions are softer in 
the chiral model, since they are coupled to constituent quarks that carry only a fraction 
of the nucleon momentum. The x dependence of the E866 data favors the virtual-pion 
model over the chiral model, suggesting that correlations between the chiral constituents 
should be taken into account. 



Recently, the flavor asymmetry of the nucelon sea was computed in the large- Ai'c limit, 
where the nucleon is described as a soliton of an effective chiral theory |l^, In this 
chiral quark-soliton model, the flavor non-singlet distribution, d{x) — u{x), appears in the 
next-to-leading order of the l/N^ expansion |ll[. The E866 d{x) — u{x) data were shown 
to be well described by this model P| . 



Instantons have been known as theoretical constructs since the seventies 45, 46, 47 



They represent non-perturbative fluctuations of the gauge fields that induce transitions 
between degenerate ground states of different topology. In the case of QCD, the collision 
between a quark and an instanton flips the helicity of the quark while creating a qq pair 
of different flavor. Thus, interaction between a u quark and an instanton results in a m 
quark of opposite helicity and either a dd or ss pair. Such a model has the possibility 
of accounting for both the flavor asymmetry and the "spin crisis" [48, 49]. However, 
the prediction at large x, d{x)/u{x) — > 4, is grossly violated by experiment (see Fig- 



ure]^). Thus, it appears that while instantons have the possibility for accounting for flavor 
and spin anomalies, the approach is not yet sufficiently developed for a direct comparison. 
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The interplay between the perturbative and non-perturbative components of the nu- 
cleon sea remains to be better determined. Since the perturbative process gives a sym- 
metric d/u while a non-perturbative process is needed to generate an asymmetric d/u sea, 
the relative importance of these two components is directly reflected in the d/u ratios. 
Thus, it would be very important to extend the Drell-Yan measurements to kinematic 
regimes beyond the current limits. 

The 50-GeV PS presents an excellent opportunity for extending the d/u measurement 
to larger x {x > 0.25). As mentioned earlier, for given values of xi and X2 the Drell-Yan 
cross section is proportional to 1/s, hence the Drell-Yan cross section at 50 GeV is roughly 
16 times greater than at 800 GeV. Figure ^ shows the expected statistical accuracy for 
(t(p -t- d)/2a{p + p) at the 50-GeV PS (see Section 3) compared with the data from E866 
and a proposed measurement using the 120 GeV proton beam at the Fermilab Main- 
Injector. A definitive measurement of the d/u over the region 0.25 < x < 0.7 could indeed 
be obtained at the 50-GeV PS. 



2.4 Polarized Drell-Yan at the 50-GeV PS 

Despite extensive work on polarized DIS, the helicity distributions of u and d sea quarks 
are still poorly known. Both the SMC and the HERMES experiments attempted 
to extract the sea-quark polarizations via semi-inclusive polarized DIS measurements, and 
the results indicate small sea-quark polarization consistent with zero. However, as pointed 
out in Ref. 
extraction. 



54 , large uncertainties are associated with certain assumptions made in the 



A direct measurement of sea-quark's polarization is clearly very important for un- 
derstanding the flavor decomposition of proton's spin. Different theoretical models make 
drastically different predictions. In particular, the meson-cloud models, which successfully 
describe the unpolarized d/u asymmetry, predict negligible amount of sea-quark polariza- 
tion |56|. Several current parametrizations |5^, ^ of polarized parton distributions 



also assume very small polarization for sea quarks. The chiral-quark soliton model, on the 
other hand, predicts substantial sea-quark polarization |11, 54]. Figure ^ shows xAu{x), 
xAd{x), and xAs{x) at Qq = 0.36 GeV^ from a recent prediction of chiral-quark soliton 
]59| . Also shown in Figure ^ are the GRSV parametrizations |Q from a global fit 



model 



to polarized DIS data. 



A very striking prediction of the chiral-quark model is the large flavor asymmetry of 
polarized sea-quark polarization. In fact, this model predicts a significantly larger values 
for Au — Ad than for d — u. This is shown in Figure |^, where x{Au — Ad) from the chiral- 
quark soliton model is compared with the x{d — u) parametrization from GRV94 . 



Polarized proton beam at the 50-GeV PS would offer an exciting opportunity for 
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probing sea-quark polarizations. The longitudinal spin asymmetry in the DY process is, 
in leading order, given by []6l| . 



Eg el[Aqa{xi)Aqa{x2) + Ag„(a:i)Ag,(x2)] 
Ea el[qaixi)qaix2) + qaixi)qaix2)] 



with Aqa = Qa ~ la ■ The superscripts refer to parton spin projections parallel (+) or an- 
tiparallel (— ) to the proton's spin projection. We have simulated the performance of the 
proposed high-mass dimuon spectrometer for measuring polarized antiquark distribution. 
Figure || shows the X2 dependence of AfJ, integrated over the spectrometer acceptance, 
for polarized sea-quark parametrizations including Gehrmann- Stirling (G-S) sets A and 
C and GRSV Leading-Order set [58|. Very small values for A^^ are predicted for the 
G-S parametrization, while the GRSV parametrization gives All ~ -0.2. The chiral- 
quark soliton model gives large positive A^}^ . In fact, the positivity requirement, namely, 
— 1 < Au{x)/u{x) < 1, is not always satisfied at the region x > 0.2 for the particular 

m. 



parametrization given by Ref. 



We have calculated the expected statistical sensitivities for a 120-day p + p measure- 
ment, assuming 75% polarization for a 5 x 10^^ per spill polarized proton beam. We 
also assume a polarized solid NH3 target similar to the one used by the SMC which 
achieves a hydrogen polarization of 75% and a dilution factor of 0.15. The target length is 
chosen to give the same gm/cm^ as for the liquid deuterium target. Figure § shows that 
the statistical accuracy of such a measurement can well test the predictions of various 
model (note that the chiral-quark soliton model predicts a large positive ^4^2^ not shown 
in this figure). A comparison of p + p with p + d will further determine Ad, which provides 
a direct test of the chiral-quark soliton model's prediction of large Au — Ad. 



In the polarized Drell-Yan process one may also measure a new structure function, 
called transversity, which is a correlation between quark momentum and its perpendicu- 
lar spin component ||63|. The transversity is not measurable in inclusive DIS 1^^. It is 
measurable, in principle, in collisions of polarized protons whose spins are aligned per- 
pendicular to the plane of dilepton detection [^. A non-zero transverse spin correlation 
in the Drell-Yan process would clearly require both quark and antiquark transversities to 
be non-zero. Polarized proton beam at the 50-GeV PS could provide unique information 
on the transversities at large x. 



2.5 Nuclear effects of Drell-Yan 



67, 



that this effect 



Following the discovery of the EMC effect, it was suggested 
is caused by the excess of virtual pions in nuclei, which significantly modify the nuclear 
parton distributions. A direct consequence of the "pion-excess" model is the nuclear en- 
hancement of antiquark sea, which can be probed via Drell-Yan experiment |^|. However, 
the subsequent Fermilab E772 experiment found no evidence for such enhancement 



[see Figure ^). The lack of an antiquark enhancement in nuclei suggests that there are 
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no more pions surrounding an average nucleon in a heavy nucleus than there are in a 
weakly bound system, deuterium. This contradicts conventional wisdom and is also at 



odds with sophisticated calculations using realistic nuclear force [70]. Unfortunately, the 
error bars for the E772 data in the region a; > 0.15 become quite large, due entirely to 
limited statistics. Furthermore, at x < 0.1 the on-set of the shadowing effect makes the 
isolation of possible pion-excess effect somewhat uncertain. 

At the 50-GeV PS, one can measure the nuclear effect over the large x region {x > 0.15) 
with high accuracy. This is illustrated in Figure where the expected statistical errors 
for a 60-day measurement of p+Ca and p+d using the proposed spectrometer (see Section 
3) are shown. One advantage of the 50-GeV measurement is that shadowing effect is no 
longer important at large x. The precise measurement at x larger than E772 could access 
would provide extremely valuable new information on the nuclear dependence of antiquark 
distributions. The anticipated sensitivity will be sufficient to observe the reduction in the 
nuclear sea distributions predicted in the rescaling models |71|. The pion-excess model. 



on the other hand, predicts a strong nuclear enhancement of Drell-Yan cross sections in 
this X region. 

2.6 Partonic energy loss in nuclei 

The subject of energy loss of fast partons propagating through hadronic matters has 
attracted considerable interest recently |]72[. The nuclear dependence of the Drell-Yan 
process provides a particularly clean way to measure the energy loss of incident quarks 
in a cold nuclear medium. Partonic energy loss would lead to a degradation of the quark 
momentum prior to annihilation, resulting in a less energetic muon pair. Therefore, one 
expects the Drell-Yan cross sections for heavier nuclear targets to drop more rapidly at 
large Xi {or xp). 



Data from E772 at 800 GeV/c were analysed by Gavin and Milana ||73| to deduce 
the initial-state quark energy loss. They ignored the shadowing effect and assumed the 
following expression for the average change in the momentum fraction: 

Axi = -KiXiA^/^. (9) 

A surprisingly large fractional energy loss (~ 0.4%/ fm) was obtained. This result 
was questioned by Brodsky and Hoyer |j7^, who argued that the time scale for gluon 



bremsstrahlung need to be taken into account. Moreover, as pointed out in Ref. [^, it is 
important to account for the shadowing effect before a reliable value of partonic energy 
loss can be extracted. Using an analogy to the photon bremsstrahlung process, Brodsky 
and Hoyer suggested an alternative expression: 

Axi ^ --A'/\ (10) 
s 

where s is the square of the nucleon-nucleon center-of-mass energy. Note that Eq. | im- 



plies a linear dependence of the energy loss on the partonic energy, while Eq. |I0| assumes 
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a constant energy loss independent of the partonic energy (note that AE is proportional 
to Axis). Based on uncertainty principle, Brodsky and Hoyer concluded that energetic 
partons should loose < 0.5 GeV/fm in nuclei. More recently, Baier et al. and Za- 
kharov [jl4| predicted 

Aa;^ ~ -1lA^/\ (11) 
s 

These authors obtained the nonintuitive result that the total energy loss is proportional 
to the square of the path length traversed. 



Very recently, the E866 nuclear-dependence data have been analysed by taking into 
consideration the shadowing effect and comparing with the three different expressions 
(Eqs. H - |n|) for energy loss []TB|. Upper limits of ^2 < 0.75 GeV^ and Kg < 0.10 GeV^ 
have been obtained. The ^2 limit corresponds to a constant energy loss rate of < 0.44 
GeV/fm, while the limit implies AE < 0.046 GeV/fm^ x L^, where L is the quark 
propagation length through the nucleus. This is very close to the lower value given by 



Baier et al. 13 for cold nuclear matter. 



A much more sensitive study of the partonic energy loss could be carried out at the 
50-GeV PS. We have simulated the effect of initial-state energy loss on the p + W Drell- 
Yan cross sections, and the results are shown in Figure |10[ Assuming a 60-day run with 
the nominal spectrometer configuration (see Section 3), the expected xi distribution for 
p + d is shown as the solid curve. The dashed, dotted, and dash-dotted curves in Fig- 
ure correspond to p + W Xi spectra assuming a partonic energy loss form of Eq. |10] 
with dE/dz of -0.1, -0.25, -0.5 GeV/fm, respectively. The ratios oi p + W over p + d, 



shown in Figure |iO], are very sensitive to the quark energy loss rate, and the expected 
statistical accuracy can easily identify an energy loss as small as 0.1 GeV/fm. The greater 
sensitivity at 50 GeV is due to the 1/s factor in Eq. |10| and Eq. |lT[ Another important 
advantage at 50 GeV is the absence of shadowing effect, and no shadowing correction is 
required. 



The Drell-Yan A-dependence data could further be used to determine whether the 
energy loss follows an L (as in Eq. |TD|) or an Li^ (as in Eq. |rTp dependence. This is 
illustrated in Figure |ll], where the solid circles correspond to {p + A)/{p + d) assuming 
an energy-loss rate of 0.25 GeV/fm using Eq. |10[ The open squares correspond to the 
situation when energy loss is described by Eq. |TT] (the value of is selected by matching 
the {p + W)/{p + d) values for both cases) . Figure |ll| shows that one can easily distinguish 
an L- from an L^-dependence even when the energy loss rate is as small as 0.25 GeV/fm. 



2.7 Quarkonium Production at 50 GeV 

Unlike the Drell-Yan process, the mechanisms for J/^ production are not well understood. 
Several quarkonium production models have been considered in the literature, including 
color-evaporation, color-singlet, and color-octet models. For simplicity, we consider the 



12 



color-evaporation model, which is capable of describing the energy-dependence and the 
shape of the differential cross sections well. However, the absolute normalization of this 
model is treated parameter. 



Figure [T^ shows the prediction of the color-evaporation model for J/^ production at 
50 GeV. The absolute normalization is obtained from an extrapolation of the global fit 
of existing J/\l' data [^. Unlike the situation at 800 GeV where the gluon-gluon fusion 



subprocess dominates ||76[, Figure ^ shows that the quark-antiquark annihilation is the 
dominant subprocess at 50 GeV. While this is reminescent of the Drell-Yan process, it 
is worth noting that quarkonium production is a hadronic process unlike the electromag- 
netic Drell-Yan process. Hence, there is no weighting factor for the q — q subprocess 
in J/\E' production. 

As indicated in Figure 0, the J/\E' production data at 50 GeV are largely sensitive to 
quark distributions and could provide information similar to Drell-Yan. This is illustrated 
in Figure ITB| which shows that the J/\E' cross section ratio for p + d over p + p is very 



sensitive to the d/u asymmetry just like the Drell-Yan process. This could be readily 
tested at the 50-GeV PS, since the J/"^ event rate is expected to be very high. 

It is also of interest to measure J/\I' production using polarized proton beam and tar- 
get. Unfortunately, the uncertainty of the production mechanism might make it difficult 
to deduce information on polarized structure functions. 



3 Experimental Apparatus 

The spectrometer considered here is designed to measure muon pairs at M^+^- > 1 GeV 
with 50 GeV proton beam. The E866 spectrometer and its daughter, a proposed P906 
spectrometer are taken as a starting point. The E866 spectrometer is shown in 

Figure n. fi~^fi~ pairs produced at the target were analyzed by a vertical-bending spec- 
trometer. The remaining proton beam was intercepted by a copper beam dump located 
inside the dipole magnet. The beam dump was followed by a set of absorbers made of 
copper, carbon, and polyethylene, which absorbed many of the pions and kaons produced 
at the target before they could decay into muon backgrounds. Trigger hodoscopes, muon 
identifiers, and tracking counters followed the magnets. The magnetic current can be 
adjusted to optimize the acceptance of a selected mass range. The spectrometer has good 
acceptance for dimuons with xp > and pt up to 3 GeV/c. 

To design a spectrometer suitable for 50 GeV proton beam, it is useful to consider 
some kinematics of the Drell-Yan process. Table |l| compares the total center-of-mass 
energy and the Lorentz factor for proton beams of 50 GeV (at the present project), 120 
GeV (at Fermilab Main Injector) and 800 GeV (at Fermilab Tevatron). 
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Table 1: The center-of-mass energy and the Lorentz factor for three beam energies 





50 GeV 


120 GeV 


800 GeV 


7/ 


9.865 GeV 
5.259 


15.12 GeV 
7.998 


38.79 GeV 
20.65 



The Lorentz factor of the nucleon-nucleon center-of-mass frame is 



El +7112 

7/ = 7^- (12) 



For fi'^ and fi emitted at 90° in the nucleon-nucleon center-of-mass frame and for Xi? ~ 
the opening angle 6 of the two muons in the laboratory frame is expressed as 



tan(^/2) = — . (13) 

7/ 

For muons emitted at 90° in the virtual-photon rest frame, the laboratory kinematics of 
the muons largely depends on M and X2 and does not depend on beam energy (or y/s). 
More specifically, 

pif = M/2. (14) 



j^-b ^ I jMb ^ I beam ^ ^ ^ (15) 



The measurable pt range of muons should be almost the same as the E866 spectrom- 
eter or slightly smaller because the mass (M^+^-) range interested is the same or slightly 
less. Thus the total magnetic rigidity (/ Bdl) of the magnets should be about 7 to 10 T-m, 
though it depends on the geometrical layout of the detectors. The P906 spectrometer has 
a magnetic rigidity of 8 T-m. 



According to Eq. |T3], the opening angle of the muons at 50 GeV is about 4 times 
larger than at 800 GeV (E866), and 1.5 times larger than at 120 GeV (P906). One idea 
to design the 50-GeV spectrometer is just to shorten the existing setup in the beam (z) 
direction with the factor of 7(50 GeV)/7(800 GeV) or 7(50 GeV)/7(120 GeV). However, 
since the maximum field strength of the P906 magnet is already near the saturation point, 
a magnet of almost the same length and wider aperture need to be considered. 



Figure |I5| shows the horizontal and vertical view of the proposed spectrometer. The 
target is assumed to be a liquid hydrogen or deuterium target 20-inch long and 3-inch 
wide. The produced charged particles are analyzed by a vertical-bending magnet, which 
is basically the same as the P906 magnet but has a wider aperture. The length of the 
magnet along the beam axis is 480.06 cm. The horizontal gap of the magnet at the exit 
is 116.84 cm and the vertical gap at the exit is 279.4 cm. The total momentum kick by 
this magnet is about 2.5 GeV/c. The incident proton beam is stopped by a copper beam 
dump, followed by a set of absorber materials. The second magnet and detectors are 
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placed after the first magnet. The momentum kick by the second magnet is 0.5 GeV/c. 
The total length of the spectrometer system from the entrance of the first magnet to the 
end of the detector system is 1474.47 cm. The muons which hit all the detectors are 
accepted as signals. 

A fast Monte-Carlo code, which takes into account the Drell-Yan cross section and the 
spectrometer configuration, has been used to estimate a statistical error for a{pd) /2a{pp) 
shown in Figure ^. In order to estimate the yields and statistical errors, the following 
assumptions have been applied: 

• The beam intensity is 1 x 10^^ protons/(3 sec). 

• The net efficiency of data acquisition is 0.5. 

• Data are taken for 60 days each for 50-cm long proton and deuteron targets. 

The performance of the spectrometer ioi p + p measurement, nuclear dependence study 
of Drell-Yan, and J/'^ production, has also been simulated and the results have been 
presented in the previous Section. 

4 Summary 

We present a broad range of physics topics which can be pursued at the 50-GeV PS 
using a dimuon spectrometer and a primary proton beam of 10^^ per spill. The expected 
sensitivities of various measurments have been simulated for a preliminary design of the 
dimuon spectrometer. The physics scope can be considerably enlarged with the addition of 
polarized proton beam and with heavy-ion beams. More detailed studies using GEANT- 
based simulation are in progress to address the issues of background and to optimize the 
design of the spectrometer. Based on our study thus far, it is clear that a rich physics 
program can be mounted using the primary proton beam at the 50-GeV PS. 
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Figure 1: Proton- induced Drell-Yan production from experiments NA3 ||25| (triangles) at 
400 GeV/c, E605 |2| (squares) at 800 GeV/c, and E772 ^ (circles) at 800 GeV/c. The 
lines are absolute (no arbitrary normalization factor) next-to-leading order calculations 
forp + d collisions at 800 GeV/c using the CTEQ4M structure functions [p8 |. 
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Figure 2: Comparison of Drell-Yan cross section data with NLO calculations using 
MRST 13 structure functions. Note that r = XiX2. The E772 [|7|, E605 m, and 
NA3 data points are shown as circles, squares, and triangles, respectively. The solid 
curves correspond to fixed-target p + d collision at 800 GeV, while the dashed curve is for 
p + d collision at 50 GeV. 
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Figure 3: The ratio of d/u in the proton as a function of x extracted from the Fermilab 
E866 P] cross section ratio. The curves are from various parton distributions. Also shown 
is the result from NA51 [0, plotted as an open square. 
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Figure 4: Comparison of the E866 |§| d—u results at = 54 GeV^/ with the predictions 
of pion-cloud and chiral models as described in the text. The data from HERMES 
are also shown. 
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Figure 5: {p + d)/{p + p) Drell-Yan ratios from E866 (open circles) are compared with 
the expected sensitivites at the 120 GeV Main Injector (sohd circles) and the 50-GeV PS 
(solid squares). 
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Figure 6: Sea quark polarizations at Q^=0.36 GeV^ predicted by chiral-quark soliton 
model are shown as solid curves. The dashed curve corresponds to the parametrization 
of x^q at g2 = 0.4GeV2 by GRSV, where Ag = An = AJ = As. 



26 



0.15 



0.1 - 



0.05 



-0.05 



-0.1 - 



-0.15 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



X 



Figure 7: x{/S.u — AJ) at = 0.36 GeV^ predicted by the chiral-quark soliton model is 
shown as the sohd curve. The GRV94 LO parametrization of x{d — u) at = 0.4 GeV^ 
is shown as the dashed curve. 
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Figure 8: Expected statistical accuracy for measuring the double-helicity asymmetry ^^2" 
in polarized p + p Drell-Yan at the 50-GeV PS for a 120-day run. The dashed, dotted, 
and dash-dotted curves correspond to calculations using polarized PDF parametrization 
of G-S (set A, set C) and GRSV, respectively. 
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Figure 9: {p + Ca)/{p + d) Drell-Yan ratios from E772 (open circles) are compared with 
the expected sensitivites at the 120 GeV Main Injector (sohd circles) and at the 50 GeV 
PS (sohd squares). 
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Figure 10: a): Solid curve is the expected p + d spectrum for a 60-day run at 50 GeV. The 
dashed, dotted, and dash-dotted curves correspond to p + W spectra assuming a partonic 
energy loss rate of 0.1, 0.25, 0.5 GeV/fm, respectively, b): Sohd circles show the expected 
statistical errors for {p + W) /{p + d) ratios in a 60-day run for p + W and p + p each. The 
solid, dashed, and dotted curves correspond to a partonic energy loss rate of 0.1, 0.25, 0.5 
GeV/fm, respectively. 30 
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Figure 11: a): Solid circles correspond to the expected {p + A)/{p + d) ratios assuming a 
partonic energy loss rate of 0.25 GeV/fm with a nuclear dependence given by Eq. TH. The 
open squares correspond to partonic energy loss given by Eq. 0. The statistical erorrs 
were calculated assuming a 60-day run for each target, b) Same as the top figure, but for 
a different Xi bin (0.7 < x < 0.8). 
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P + D ^ J/t + X at 50 GeV 
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Figure 12: Calculation of the p + d J/\& + x cross sections at 50 GeV using the 
color-evaporation model. The contributions from the gluon-gluon fusion and the quark- 
antiquark annihilation subprocesses are also shown. 
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J/t Production at 50 GeV 
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Figure 13: Calculations of the p + d —>■ J/\& over p + p —> J/\l' ratios at 50 GeV using 
the color-evaporation model. The d/M- symmetric structure functions DO 1.1 and the d/u 
asymmetric structure functions (MRST and CTEQ4M) have been used in these calcula- 
tions. 
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Figure 15: A schematic view of the prototype spectrometer. The top is the horizontal 
view and the bottom is the vertical view. 
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